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Under externally applied electric fields, lipid membranes tend to permeate and change their 
electrical resistance by the combined processes of pore creation and pore evolution (expansion or 
contraction). This study is focused on the pore creation process, represented by an empirical 
expression currently used in the electroporation (EP) models, for which an alternative theoretically 
based expression was provided. The choice of this expression was motivated by the role the DLVO's 
(disjoining) pressures may play in the process of EP. The electrostatic energy effects on each sides of 
a lipid membrane were evaluated in terms of the electrostatic component of the disjoining pressure. 
Thus the pore creation energy considerations in the current EP models, associated with the necessity 
of an idealized non conducting circular pre-pore were avoided. As a result, a new expression for 
the onset of the electroporation was proposed. It was found that this new theoretically determined 
expression is in good agreement with the one estimated from an experiment that specifically targeted 
the pore nucleation by eliminating the pore evolution process. Furthermore, the proposed expression 
is dependent on the electrolyte and the lipid properties and should provide better predictability than 
the currently used pore creation EP model, which is only temperature dependent. 

PACS numbers: 



I. INTRODUCTION 

The application of electric fields to cells with intention 
to introduce various molecules into them is widely used 
physical method. This technique of permeating the cel- 
lular membrane, called electroporation (EP), has been 
routinely performed for many decades. By applying rel- 
atively strong electric pulses, the fluid 5 nm thick bilipid 
membrane that envelops the cell and the organelles, can 
be perforated by nanometer size pores that are long-lived 
(ms) compared to the time it takes to create them (ns). 
It is still puzzling why the resealing period is at least six 
orders of magnitude longer. Most likely this is a reflec- 
tion of the difference in the physics involved in these two 
processes. These electro pores are too small to be opti- 
cally observed. The existence of these perforations in the 
membrane is usually confirmed by measuring the abrupt 
change in membrane's electrical resistance, which can be 
attributed to a combination of two different processes: 
nucleation of new pores and expansion or contraction of 
the existing ones (pore evolution). Experimentally, by 
measuring the membrane's electrical resistance, it is very 
hard to differentiate between these two processes. 

Although the main goal is to understand the electro- 
poration of cells, in order to understand the EP process 
better, much experimental and theoretical work has been 
done on planar bilayer membrane. The current theoreti- 
cal understanding of the electroporation process is built 
on the expressions originally used in the work on thin 
films. This is natural considering the composition simi- 
larity between thin films and planar bilayer membranes. 
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They both consist of a two-dimensional fluid layer of lipid 
molecules interacting with surrounding aqueous solution. 

II. MODELING THE TWO ASPECTS OF THE 
EP PROCESS - PORE NUCLEATION AND PORE 
EVOLUTION 

The basic theoretical description of hole nucleation in 
thin films goes back to 1940's in the work of Kramers 
and Zeldovich and starts with considering the ex- 
pression for reversible work of formation of a pre-pore, a 
macroscopic hole having a circular shape. The main con- 
tributors to this mechanical energy expression are the 
line tension of this circular hole and the surface tension 
of the membrane water interface. These kinds of energy 
considerations are also used in EP theory and modeling 
to describe the two main effects an applied electric field 
would have on a membrane, namely creation of new pores 
and expansion or contraction of the existing ones. 

This study does not address the process of pore evo- 
lution (expansion/contraction) in the EP models. It is 
focused only on the pore nucleation aspects of the EP 
theory and modeling. 

A. Currently used pore creation rate for the 
process of pore nucleation in the EP models 

In the current EP models the process of pore nucleation 
in membranes under applied electric field is described by 
a pore creation rate expression that besides the mechani- 
cal energy it also includes the electrostatic energy effects. 
Extending Litster's 1975 pore fluctuation model ^2], and 
Chizmadzhev 1979 work on mechanisms of defect origina- 
tion in membranes [3!l,in 1981 the following pore creation 
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rate expression was introduced by Weaver and Mintzer 
i 

dt 

with Aq — i^oAj^dia- Here N is the number of created 
pores, 1^0 is referred to as attempt rate density per vol- 
ume, Am and d,n are membrane's area and thickness, kT 
the Boltzmann's constant and the temperature, Aip^ is 
the transmembrane voltage, constant /3 in units of ca- 
pacitance is associated with the electrostatic energy ef- 
fects, and dc is the energy barrier to be overcome for 
a single pore to be formed. Thus Sc is the one pore 
formation free energy required for cooperative topolog- 
ical rearrangement of many lipid molecules from planar 
two-dimensional to most likely traversable wormhole-like 
three-dimensional structure. The actual principles of 
membrane defect formations are not known yet, but it 
is generally accepted that the process is cooperative and 
stochastic, probably involving not more than « 2x50 
lipid molecules 0]. 

A similar version of the pore creation rate was sug- 
gested in 1988 by Glaser et al. with Aq = vA^/Ahg, 
where A^g is the headgroup area of a lipid molecule 
and v is the frequency of lateral fluctuations of the lipid 
molecules. This choice of parameters in Aq is based on 
the authors of Q believe that the rate at which the pore 
can open is mainly determined by lipid headgroup fluc- 
tuations. 



B. Currently used pre factor l3 in the EP models 

Expression ([T]) includes the electrostatic energy effects 
associated with the transmembrane voltage, A^pm, which 
tend to decrease the stability of the membrane against 
thermal fluctuations. In current EP models the way 
the electrostatic energy modifies the effective mechani- 
cal pore formation energy is by considering the pore as 
a circular capacitor. The pre factor /3 (capacitance) con- 
tains the difference in energy when this capacitor is field 
with water instead of lipid membrane 

a _ 7r^» (^w ~ £m) /„x 

^" 2dm ' ^' 

where is some minimum pre-pore radius, and £w and 
Em are the corresponding dielectric constants of water 
and lipids. 

During the last two decades, most of the EP com- 
puter simulation models 0, H, H, [13, [HI [Hi use a general 
pore creation rate given by ([T]), with choice of parame- 
ters a(^'o, Aya, Sc) and P that provide a reasonable exper- 
imental predictability of the EP models. In some of these 
models [13, [Hi, it is taken that /3 = l/(Ai/)ep)^, where 
Aipep specifies the onset of the electroporation and is re- 
ferred to as the minimal transmembrane voltage that can 
cause a pore to form. 



With time, it is been always a theoretical challenge to 
accommodate the experiments demonstrating that lower 
and lower voltages can cause electroporation. Very often 
in the EP models this has been achieved by adjusting 
the value of the pre factor (3. Since /3 appears in the 
exponent, it strongly influences the situation that ([T]) is 
attempting to model. As a result, some degree of ambigu- 
ity regarding the value of /3 has been present in modeling 
the EP process. Over the past two decades the different 
values of /3 used in the EP models have ranged from 4.70 
kTY-^ to 62.50 kTY'^ (at temperature T = 293 K and 
k = 1.38 X 10-23 jK-i). 

Chronologically, in 1991 [7], by assuming minimum 
pore radius — 1 nm and membrane thickness d^ = 2.8 
nm the authors have obtained /3 — 95.70 kTY^^ but used 
a more reahstic /3 = 4.70 kTY^'^ in their EP model. By 
1996 ^, the authors considered more standardized val- 
ues of = 0.4 nm and dm = 5 nm and obtained (3 = 8.59 
kTY-"^. In 1998 the authors have assumed (3 = 62.50 
kTY^^ without providing information about the param- 
eters used to estimate this value. The authors of [l3| in 
1999 have proposed that pores are being created at all 
radii, with decreasing probability for the larger radii, and 
they suggested a value /? = 13.12 kTY~^. Similarly in 
their 1999 paper [l3|, the authors use (3 = 15.02 kTY~'^. 
But they also assume = 0.76 nm, which would give 
(3 = 30.95 The authors of ^ also take that the 

pores are being created at all the radii in a decreasing 
manner, from their minimum pore radius — 0.8 nm, 
for which f3 = 34.34 kTY-^. 

An increase in the applied voltage will result in in- 
creased membrane conductance, which can be attributed 
to a combination of two different processes. The process 
of creation of new pores in the membrane and the pro- 
cess of expansion of the existing pores. When both of 
these processes are present, experimentally it is impossi- 
ble to determine which part of the increased conductance 
is due to the newly created pores and which part is due 
to expansion of the existing ones. Thus relevant to the 
pore creation rate would be an experiment that will avoid 
the process of pore expansion and would only reveal the 
fundamentals of pore formation. 

Exactly this kind of experiment, essential for the pore 
creation rate was performed in 2001 (isj . in which a nu- 
cleation of a single pore was experimentally observed. 
This experiment was done by applying not an electric 
pulse but a constant low voltage across an A^ = 3 ^m 
X 3 /im artificial bilipid membrane patch, until a single 
pore appeared and resealed. Surprisingly, even voltages 
as low as 250 mV could cause an appearance of a sin- 
gle pore. Based on this experimental data, crucial for 
the pore creation rate, the EP model parameters were 
roughly estimated in to be a = 0.01 and /3 = 20 
kTY~^ thus giving to a and /? in expression ([1]) some 
experimental justification. 
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C. Critique of the currently used pre factor P in 
the EP models 

Further in this paper, the vahdity of the theoretical 
assuniptions, based on a pre-pore existence (i.e. circular 
capacitor), contained in expression ^ are criticized and 
an alternative expressions for the electrostatic energy ef- 
fects on pore nucleation in ^ is suggested. 

The existing theoretical approach for the pore creation 
rate is based on the reversible work, assuming the exis- 
tence of a minimum size pore in the membrane. This 
approach does not consider how the changes in the elec- 
trolyte due to the applied voltage led to the nucleation 
of this minimum size pore. As such, all of the so far 
used EP models do not take into account the electrolyte 
or lipid properties. Thus they give the same prediction 
for different types of electrolytes and counterion concen- 
trations. It is hard to imagine that solution properties 
would not play any role in the EP process. This property 
independence hinders the predictability of the EP models 
and requires recurring adjustment of the model constants 
to fit the experiments. Since constants don't change un- 
less they are parameters, for the predicting purposes of 
the EP models, it is very important to use theoretically 
justified parameters in their pore creation expressions. 

The current theoretical understanding of the expres- 
sion (2)) relies on the well-known fact that the lipid mem- 
brane behaves essentially as a local circular capacitor, 
upon the application of an electric field. Expression ([2]) 
is certainly valid for large holes in the membrane but ge- 
ometrically it cannot be circular for small ones. In this 
context all current pore nucleation models assume pre 
existence of a circular pore with a cross section TrrJ, of 
minimum pore radius r.^, « 0.4 — 1 nm. Geometrically, 
this can hardly be correct, considering that the minimum 
radius pore is made of only few (5-6) lipid hcadgroups of 
comparable radius ~ 0.8 nm (area Ahg « 0.6 nm^ [14|). 
thus making 7rr^ in ((21) a very crude approximation to a 
circular capacitor even without including the permanent 
motion of the lipid molecules. In addition, the dielectric 
properties appearing in expression ([2]) most likely differ 
from the bulk phase, taken — 80eo and £,„ = 2eo 

by all the EP models. 

In light of all this, a possible resolution to some of 
these issues related to the electrostatic energy effects in 
the pore creation rate is further suggested. 



III. APPLICATION OF DLVO ON A LIPID 
MEMBRANE 

For a number of decades the DLVO theory, started 
by Derjaguin and Landau (1941), and complemented by 
Verwey and Overbeek (1948), has been the well accepted 
theory that describes the interaction of membranes with 
the surrounding electrolyte, giving an expression for the 
(disjoining) pressure on the membrane in terms of the 
counterion properties. Thus it seems natural to apply 



DLVO to the pore creation rate used in EP models. 

A. The electric potential near a lipid membrane 

A membrane in electrolyte solution is screened from 
the environment by the electrolyte counterions that are 
distributed against it and they balance membrane's sur- 
face charge [ij, [111 • The counterion concentration in the 
vicinity of the membrane is up to two orders of magni- 
tude higher than in the bulk. This causes a side electric 
potential ips to be formed near the membrane within few 
Debye lengths (~ 1 nm) from it, as illustrated in Fig. 
[T] The magnitude of the side potential on the left t/'SL 
and the right side ipsR is determined by the membrane's 
surface charge density a, pH, as well as the bulk ion con- 
centration no in the solution on each side. According to 
Grahame Equation [3, UBIj "03 oc (jl^/n^. For instance 
if the bulk concentration is increased, the side potential 
will get smaller. This will not result in appearance of 
transmembrane voltage A?/;iu across the membrane. A 
typical value of the side potential is -67 mV, as given 
in reference [3], estimated for a 100 mM aqueous solu- 
tion of NaCl, against a typical lipid membrane surface 
charge of —0.2 Cm^^. There is an additional potential 
drop of ~ 100 mV in the Stern layer near the membrane 

Regardless of the membrane's surface charge or the 
bulk ion concentration, when no external electric field 
is applied, the electric potentials in the bulk on the left 
and on the right side are equal [3) EBl- This is illus- 
trated by the dashed line in Fig. [T] An applied electric 
field redistributes the charges and the potentials near the 
membrane and creates a transmembrane voltage (Ai/jm). 
Important to note [iJi EBl is that the externally applied 
electric field has virtually no effect on the side poten- 
tials "i/jsL and ■0SR- This indicates that the applied elec- 
tric field across a membrane manifests itself only by the 
transmembrane voltage [3, [lB| • 

B. The electrostatic component of the disjoining 

pressure 

In DLVO it is customary to separate the various con- 
tributions of the disjoining pressure into different compo- 
nents: electrostatic, dispersion, steric, adhesive, etc. In 
this paper only the electrostatic He component of the dis- 
joining pressure was considered. The expression for lie 
on one side of a single membrane exposed to electrolyte 
of charge density p, has the following form [3j [13 

n,^- J piij)d^. (3) 



The sign indicates that the force pushing the counteri- 
ons in the electrolyte is directed towards the membrane. 
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The integration by the potential tp is done at a particu- 
lar location x away from the membrane, as the potential 
^{x) at that location changes from zero to some new 
value "01. An observation made by Derjaguin [l|, [l^ in 
late 1930's illustrates that the disjoining pressure IIo, af- 
ter brief equilibration of the ions, is constant throughout 
the electrolyte (i.e. independent of the location x). Thus 
when an electric field is applied across the membrane, He 
can conveniently be estimated by the change of -0 in the 
bulk. 

In a symmetric situation, when same bulk concentra- 
tion riQ is present on each side of the membrane, and 
no external electric field is applied, the bulk electric po- 
tential is zero as illustrated by the dashed line in Fig. 
[TJ In this situation the membrane does not experience a 
difference in He on each side. 

If an external electric field is applied, it breaks the sym- 
metry in the electric potential around the membrane as 
illustrated by the full line in Fig. [TJ The potential in the 
bulk on the right increases from zero to the value of the 
transmembrane voltage A^/im. According to the defini- 
tion of Ho, expression ([3]), the membrane will experience 
a pressure Alio from the right side. 

The actual expression for the electrostatic component 
of the disjoining pressure He, for a Boltzmann charge 
density distribution p (ip) in expression ^ was derived 
in 1948 by Verwey and Overbeek. This expression ap- 
plied on one side of a single membrane exposed to 1:1 
electrolyte has the following form 



iDrane exposea 



An, 
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1 / zey 
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g2 

-no — AV'm',(4) 



where expression (|3]) was integrated over the sum of the 
anion and the cation charge density distribution (valence 
z = ±1). 

The asymmetry on each side of the membrane can be 
realized in many circumstances. For example by combi- 
nation of any of these conditions: different counterions 
(valence), bulk concentrations no, different lipid compo- 
sition (surface charge) or by applying an external electric 
field. A detailed asymmetric analysis is beyond the scope 
of this work, even though some of these asymmetric con- 
ditions may be biologically relevant. In addition there 
is no experimental data on appearance of a single pore 
under asymmetric conditions. Thus further in this paper 
only symmetric conditions were analyzed, namely same 
bulk concentration and lipid composition on each side 
of the membrane, and a comparison was made with the 
symmetric experimental set up in [l^, vital for the pore 
creation rate, since it led to appearance of a single pore. 
It also should be mentioned that the pore creation rate 
currently used in the EP models ignores any asymmetry 
since it is independent of any electrolyte properties. 



IV. DLVO BASED EXPRESSION FOR THE 
PORE CREATION RATE AND p IN THE EP 
MODELS 



Motivated by these DLVO considerations and [igl, in 
this paper, instead of the current pore creation rate given 
by ll]), the following DLVO based pore creation rate is 
proposed 



dN 
lit 



(5) 



where AH = IIl — IIr, as Fig. [T] illustrates, is the dif- 
ference in the disjoining pressure on each sides of the 
membrane. In the spirit of mean field theory, namely re- 
ducing the multi-body interaction to effective one-body 
interaction, this disjoining pressure (energy density) was 
multiplied by an effective volume aAhg of the electrolyte, 
that contributed to formation of a single pore. This is 
the volume in which a lipid hcadgroup of area Ahg feels 
the disjoining pressure from the collisions with the coun- 
terions of radius a. 

In general, the AH term in expression ([S]) contains the 
difference of all of the disjoining pressure components. 
This paper was focused only on Alio and its pre factor 
(3. Since AH is additive, the contribution of the other 
pressure components will affect only the pre factor a. 
The influence of the other disjoining pressure components 
on a in ([S]) will be analyzed in a future work. 

The general pore creation rate ([T]) was compared with 
the DLVO based pore creation rate (O , considering only 
the electrostatic component of the disjoining pressure 
(An = -Alio). Using the expression ([4]) for Alle the 
value of the DVLO based pre factor (3 was determined to 
be 



(6) 



Unlike the previously used expression ([2]) for /3 which 
requires existence of a circular pre-pore and is indepen- 
dent of the electrolyte properties, the pore creation rate 
with the new pre factor ^ is dependent on the elec- 
trolyte and lipid properties such as counterion size a, 
bulk concentration no and lipid headgroup size ^hg- In 
this context, the value of = 20 kTV~^ that was esti- 
mated in [1] is specific to the experimental setup of p^ . 
namely 100 mM bulk concentration of K+ counterions. 
An estimate of the DLVO based expression ^ , at room 
temperature T = 293 K and with counterion radius of 
aK+ = 3.3 A [11], gave a value of /? = 18.67 /cTV-^. 
This value, based entirely on theoretical grounds, is very 
close to the round, experimentally estimated one in [5|, 
bearing in mind that in /3 was estimated in conjunc- 
tion with a as a fit to the experimental data. 
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V. CONCLUSION 

In conclusion, based on DLVO theoretical consider- 
ation, an alternative pore creation rate expression was 
proposed for use in the EP models. The energy effects 
in the proposed expression are due to the difference in 
the disjoining pressure AH on each sides of the mem- 
brane. More specifically, the electrostatic energy effects 
on a lipid membrane were considered in terms of the elec- 
trostatic component of the disjoining pressure Alio, and 
a new expression for (3 was suggested that avoids the ne- 



cessity of an idealized circular pre-pore. The value of 
this new (3, founded entirely on theoretical arguments, 
was established to be in a good agreement with the esti- 
mated experimental value in ^ . Since it is dependent on 
the electrolyte and the lipid properties, this DLVO based 
j3 should provide better predictability than the current 
EP models that regard /9 as a constant. The role that 
other components of the disjoining pressure may play in 
the EP theory and how they may affect a in ([T]) will be 
considered in a future work. 
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FIG. 1: Profile of the electric potentials ^p{x) near a membrane. In the absence of an external electric field the side potential 
is represented by the dashed line on the right side. If an external electric field is applied, it will cause an appearance of a 
transmembrane voltage Aipm across the membrane, but it will have virtually no effect 1^ on the side potentials, i/jsl and ipSR, 
that depend on the lipid composition (surface charge) and on the electrolyte's bulk concentration on each side of the membrane 
0, EE]- The potential jump Aipm in the bulk, according to the definition of the electrostatic component of the disjoining 
pressure, expression ([3]), will create a pressure AH on the membrane that is given by expression Q. 



